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Abstract. - Sufficiently dimerized quantum antiferromagnets display elementary S = 1 excita- 
tions, triplon quasiparticles, protected by a gap at low energies. At higher energies, the triplons 
may decay into two or more triplons. A strong enough magnetic field induces Bose-Einstein 
condensation of triplons. For both phenomena the compound IPA-CuCls is an excellent model 
system. Nevertheless no quantitative model was determined so far despite numerous studies. Re- 
cent theoretical progress allows us to analyse data of inelastic neutron scattering (INS) and of 
magnetic susceptibility to determine the four magnetic couplings J\ ~ —2.3 meV, J2 ~ 1.2 meV, 
J3 ~ 2.9 meV and J4 ~ —0.3 meV. These couplings determine IPA-CUCI3 as system of coupled 
asymmetric S = 1/2 Heisenberg ladders quantitatively. The magnetic field dependence of the 
lowest modes in the condensed phase as well as the temperature dependence of the gap without 
magnetic field corroborate this microscopic model. 
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Low-dimensional antiferromagnetic quantum spin sys- 
tems display various fascinating properties, e.g., spin- 
Peierls transition [U[2], appearance of a Haldane gap 
for integer spins [SIS], high-temperature superconductiv- 
ity upon doping [5], and the Bose-Einstein condensation 
(BEC) in spin-dimer systems jBHS]? where the latter one is 
characterized by a phase transition from a non-magnetic 
phase to a long-range antiferromagnetically ordered gap- 
less phase at a critical magnetic field H&. 
' Another fascinating phenomenon recently observed in 
low-dimensional antiferromagnets is the decay of their el- 
ementary 5 = 1 excitations, triplons [10], at higher ener- 
gies so that the triplons exist only in a restricted part of 
the Brillouin zone JTTJ[12] . Theoretically as well, there is 
rising interest in the understanding and quantitative de- 
scription of this phenomenon for gapped triplons [T3UT6] 
as well as for gapless magnons [T7HT9] . 

The description of quasiparticle decay faces an intrin- 
sic difficulty. The merging of the long-lived, infinitely 
sharp elementary triplon with a multitriplon continuum 
requires to describe the resulting resonance and its edges 
precisely. This is still a challenge for numerical approaches 
such as exact diagonalization or dynamic density-matrix 
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renormalization [20]. Diagrammatic approaches are able 
to capture the qualitative features but may encounter dif- 
ficulties in the quantitative description in the regime of 
strong merging where the sharp mode dissolves completely 
in the continuum because this is a strong coupling phe- 
nomenon [13l[T3]. Unitary transformations also face diffi- 
culties when modes of finite life-time occur [T6] . 

A crucial step in the understanding of both phenomena 
is to identify a suitable experimental system. The best 
studied candidate for the BEC in coupled spin-dimer sys- 
tems is TICUCI3. Unfortunately, recent research suggests 
that the high field spectrum remains gapped [HI [22] in 
contrast to what is expected from a phase where a con- 
tinuous symmetry is broken. This suggests the existence 
of anisotropics. A promising alternative for a BEC in a 
spin-dimer system is (CHa^CHNHaCuCla (isopropylam- 
monium trichlorocuprate(II), short: IPA-CUCI3) where in- 
elastic neutron scattering (INS) provides evidence for an 
almost exact realization of a BEC [91123]. 

A suitable experimental system to study triplon decay 
in detail is searched for. The two-dimensional (2D) PHCC 
[TTJ[23] i s a candidate, but it involves eight different cou- 
plings so that a quantitative characterization is impossi- 
ble to date. Due to its quasi one-dimensional (quasi ID) 
structure, IPA-CUCI3 is again a more promising candidate. 
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This compound seems to realize the theoretically proposed 
situation for BEC in coupled spin ladders [8]. 

But in spite of many years of intensive studies [9|fT2| 
[251429] no quantitative microscopic model for IPA-CUCI3 
is established. The present work aims at filling this gap. 
Theoretically, our study is based on continuous unitary 
transformations (CUTs) of models with quasiparticle de- 
cay [16] and on high temperature series expansions for 
asymmetric spin ladders which are topologically equiva- 
lent to dimerized and frustrated spin chains [30] . The 
experimental input used in INS data [12] and magnetic 
susceptibility x(^) data [26]. We will illustrate why it is 
intrinsically difficult to determine the microscopic model. 

Finally, we will compute the temperature and the mag- 
netic field dependence of the lowest magnetic modes as 
well as the upper critical magnetic field H C 2 , which induces 
full polarization. They all agree very well with experimen- 
tal data [9 ,23 ,31-33] which supports the advocated model. 
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Fig. 1: (Color online) Sketch of IPA-CUCI3. Circles indicate Cu 
ions with S = 1/2. The couplings Ji and J4 are ferromagnetic 
(Ji, J a < 0) while J2 and J3 are antiferromagnetic (J2, J3 > 0). 
Two spins linked by J 3 form a dimer. 

Since the characterization of IPA-CUCI3 by Roberts et 
al [25] various spin models were discussed. Manaka et 
al pointed out that the magnetic susceptibility of IPA- 
CUCI3 can be explained by a ferro-antiferromagnetically 
alternating Heisenberg 5 = 1/2 chain with ferromagnetic 
coupling twice as large as the antiferromagnetic coupling 
[26] . According to Hida [3l] the magnetic ground state is 
thus given by a gapped Haldane state [3] . 

The dispersions measured by INS [12] and the crystal 
structure of IPA-CUCI3 indicates that the system is quasi- 
2D. It is described by weakly coupled asymmetric spin 
5=1/2 Heisenberg ladders, see Fig. [T] with 

H =H 1D + H ± (la) 

r,s r,s 

+ J2 ^ (Sl 5r , s Si ?r+ i ?s + S2,r,sS2,r+l,s) (lb) 

r,s 

H± =J4 Si >r>s S2,r+l,a+l (lc) 

r,s 



with two ferromagnetic couplings Ji, J4 < and two 
antiferromagnetic couplings J2, J3 > 0. The dominant 
dimer coupling is J3 so that we use the ratios x = J2/J3, 
y = J1/J3 and z = J4/J3. Let us first consider the ladders 
as isolated because the interladder coupling is small. The 
standard view of these ladders takes the J3 bonds to form 
the rungs of the ladder. Then J\ is a diagonal bond. 

The key element of this model is the asymmetry of the 
spin ladders controlled by J\. On the one hand, the pres- 
ence of Ji spoils the reflection symmetry about the center 
line of the ladder between the legs. This symmetry would 
imply a conserved parity such that the triplons on the 
dimers could be changed only by an even number [35]l36] 
so that no decay of a triplon into a pair of triplons could 
occur. Hence the very presence of J\ opens an important 
decay channel for quasiparticle decay. 

On the other hand, the two bonds J2 and J\ represent 
the coupling of adjacent dimers. Both contribute to the 
hopping of the triplons which is given in leading order by 
2 J 2 — J\ [37] while the interaction of adjacent triplons is 
proportional to 2J2 + J\. With information only on the 
dispersion [12] it is impossible to determine J\ and J2 
separately. Hence, the same feature that induces the in- 
teresting quasiparticle decay makes it particularly difficult 
to establish a microscopic model. 

The BEC occurring in TICUCI3 was successfully de- 
scribed by the bond-operator approach [38,39 . But this 
approach to spin-dimer systems is quantitatively reliable 
only as long as the interdimer couplings Jmter are signif- 
icantly smaller than the dimer coupling Jdimer- | Winter | < 
^dimer/2 [40 j . This limit requires \Ji\ < J3/2 for i E 
{1,2,4} for IPA-CuCl 3 which does not hold |2B]. We will 
see below that |Ji| w J3 provides very good fits. 

Thus we apply self-similar CUTs (sCUTs) to isolated 
ladders [4TH44] , modified to cope with decaying quasipar- 
ticles [16. . We use an infinitesimal generator which decou- 
ples the subspaces with zero or one triplons from the re- 
maining Hilbert space. We can still decouple the 1-triplon 
subspace from the 2-triplon subspace for the isolated lad- 
der. The proliferating flow equations are truncated if the 
range of the corresponding process exceeds certain maxi- 
mum extensions in real space 0. Thereby, the ladders are 
mapped to an effective model 

#lD,eff = ^(WlhSohJ (2) 

h,l;a 

in terms of triplon creation h l and annihilation opera- 
tors t a h l in momentum space, where h is the wave vector 
component along the ladders, I the one perpendicular to 

1 The truncation scheme used for the Hamiltonian is 
(cfe, ^3, . . . , ds) = (10,8,8,5,5,3,3) and for the observables 
(di, d,2, . . . , do) = (10,10,8,8,6,6), where dj is the maximum 
extension for a process with j creation and annihilation operators. 
Additionally, we keep only terms that create or annihilate at most 
N = 4 triplons in the Hamiltonian and N = 3 triplons in the 
observables, see also Ref. [T6ll43| . 
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them, and a G {x, y, z} the spin polarization. These op- 
erators are the Fourier transforms of the bond operators 
[13I16] defined on the dimers in Fig. [TJ 

The dispersion u)o(h) depends only on h because the 
CUT is applied to the isolated ladders which still have to 
be coupled. This coupling is achieved in leading order fol- 
lowing the approach in Refs. [37JHH]. The spin component 
Sf r s is taken as observable and transformed into the new 
basis by the CUT. Then it reads 



a,r+5,s 



) + -.., (3) 



where the dots stand for normal-ordered higher terms in 
the real space triplon operators t\ xrs (t a ^ s ). Knowing 
^eff i r s allows us in a second step to write down the effec- 
tive interladder coupling i/mt,eff in real space 



#int,eff - Ja ^2 5Z ai ' 5a2 ' 5/ [^^^(^,r+l + (^-5), S +l 



r,s;a 5,5' 



,)+H.c.]. 



(4) 



*,r+l+(<S'-<S),s+l' 

This neglects trilinear and higher contributions. The 
Fourier transform of i?int,eff leads to H e R = i^iD,eff + 
^int,eff amenable to a Bogoliubov diagonalization yielding 



H e V=^Uj(h,l)bl hl b ahl 

h,l;a 



u(h,l) =\Jul(h) + 4cu (h)\(h, I) 

X(h, I) = — J4 ^ CLi^aij' cos (27r [h (5 + 5 f 



1) 



(5a) 

(5b) 
I]) 



5,5' 



with bosonic operators b J a h l (b a ^,i). I n the Bogoliubov 
diagonalization the hardcore property of the bosons is ne- 
glected. However this does not concern the large intralad- 
der couplings, but only the small interladder couplings so 
that the approach is still very accurate 05]. The disper- 
sion uj(h,l) makes a direct comparison with INS results 
possible. 

Table 1: Parameters for IPA-C11CI3 compatible with INS Q2] 



J3 [meV] 


x J2/J3 


y J1/J3 


z J4/J3 


3.743 


0.133 


-2.0 


-0.076 


3.288 


0.268 


-1.4 


-0.088 


3.158 


0.317 


-1.2 


-0.092 


3.038 


0.369 


-1.0 


-0.096 


2.929 


0.424 


-0.8 


-0.100 


2.830 


0.480 


-0.6 


-0.103 



To determine the microscopic parameters we fix the 
value y = J1/J3 and fit x = J2/J3, z = J4/J3, and the 
energy scale J3 to reproduce the experimental result (Eq. 
(2) in Ref. [12]) 

uo{h,l) 2 = a 2 cos 2 (irh) + [A 2 + 4b 2 sm 2 (irl)} sm 2 (irh) 
+c 2 sin 2 (27r/0 (6) 



with a = 4.08(9)meV, A = 1.17(l)meV, b = 0.67(l)meV 
and c = 2.15(9)meV. Thus, we obtain the triples (x,y,z) 
in Tab. [1] They all essentially imply the same dispersion, 
see Fig. [2j Hence, on the basis of the the INS data, one 
cannot decide which of the triples applies to IPA-CUCI3. 

The quasiparticle decay occurs where the dispersion en- 
ters the 2-triplon continuum. It does not prevent to use 
the CUT for the isolated ladder since the realistic param- 
eters turn out to be such that the triplons do not decay 
without the interladder coupling. A quantitative descrip- 
tion of the decay is subject of ongoing research. 
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Fig. 2: (Color online) Circles are INS data [12]. (a) Dispersions 
cj(h, 0) for xyz triples in Tab. [TJ The quasiparticle decay oc- 
curs where the dispersion enters the 2-triplon continuum, (b) 
Dispersion cj(0.5, /); all triples lead to coinciding curves. 

In complement to the INS we use the temperature de- 
pendence of the magnetic susceptibility x(T) [26]. Start- 
ing from the spin isotropic Hamiltonian (pQ) the susceptibil- 
ities in different spatial direction have to be the same up to 
scaling proportional to the squares of the Lande ^-factors. 
This means that Xa • Xb • Xc equals g\\ g\\ g 2 ^ where A, 
B, C indicate the directions normal to the corresponding 
surfaces of the crystal [26J. Fig. [3^ displays that the three 
susceptibilities can be scaled to coincide for gx = 2.08, 
#b = 2.06, and gc = 2.25 within about 3%. This choice 
of ^-factors fulfills the experimental constraints [26] [27] 
g A ,g B e [2.06,2.11] and g c = 2.25 - 2.26 best. We con- 
clude that an spin isotropic Hamiltonian such as (pp) pro- 
vides a very good description, although anisotropics, e.g., 
Dzyaloshinskii-Moriya terms, can be present with a rela- 
tive size of a few percent. This agrees with findings from 
electron paramagnetic resonance [27] . 

Theoretically, we use the high temperature series ex- 
pansion for the isolated asymmetric ladder [30] provid- 
ing series in f3 = 1/T up to order j3 n+1 with n = 10 
denoted by xid- The 2D series X2D obeys the relation 
X2D = Xid + ^ 4 m interladder mean-field approximation, 
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Fig. 3: (Color online) Upper panel: Deviations of the experi- 
mental magnetic susceptibilities [26] xb and xc in B and C di- 
rection relative to xa for gx — 2.08, #b = 2.06, and gc — 2.25, 
indicating anisotropics. Lower panel: Comparison of Xa(T) 
for various values gx with theoretical results obtained by Dlog- 
Pade approximated high temperature series expansions for the 
xyz triples from Tab. [T] 



i.e., in leading order in J4. We use standard Dlog-Pade 
approximation [50] to deduce the full x(T) from X2D and 
from the asymptotic behavior X2d(/?) /3° exp (— A/3) for 
<C A. The result! is plotted in Fig. [3] and compared 
to Xm measured in [emu/g] and converted according to 
X(T) = m mo ik B (gmy 2 N^XmiT). Here m mo i is the mo- 
lar mass of IPA-CUCI3, &b the Boltzmann constant, jub the 
Bohr magneton and TVa the Avogadro constant. 

Fig. [3}d illustrates that theory and experiment agree in- 
deed best for gx = 2.08 and the triple of y = —0.8. As 
an asset, we stress that even without the value of #a, the 
position and the shape of the maximum of xCO fits best 
for the triple of y = —0.8 and one can deduce deduce 
that the #A-factor is around 2.08. As a caveat, we stress 
the very weak dependence of x(^) on V m a triple tuned 
to the INS data. By assuming gx = 2.08 ± 0.01 we es- 
timate the error of our analysis to be x = 0.42 ± 0.06, 
y = -0.8 ± 0.2 and z = -0.100 ± 0.004 implying J x = 
-2.3 ± 0.6 meV, J 2 = 1.2 ± 0.2 meV, J 3 = 2.9 ± 0.1 meV 
and J4 = —0.292 ± 0.001 meV. These values establish the 
microscopic model for IPA-CUCI3. We highlight that the 
ferromagnetic coupling J\ does not dominate over the an- 
tiferromagnetic coupling J3 because \y\ £ 1, in contrast to 

2 All theory curves rely on the [7,4] Dlog-Pade approximant in 
u = /3/(l + /3). Data from other Dlog-Pade approximants, e.g., [9,2], 
agrees within line width except at very low temperatures. 



the previous purely ID analysis [26] . 

The derived microscopic model successfully passes three 
checks: The BEC is well-described, the upper critical field 
H C 2 agrees to experiment and the temperature dependence 
of the spin gap matches recent data. 

First,we follow Refs. [38|l39|f5Tj to describe the BEC and 
perform the local transformation 



iQor 



S T ) 



\sT)=u\s r )+ve^(f\U, r )+g\t_, r )) 

|?+,r) =u(f |t+, r ) + g |£_ r )) - 
|^0, r) = |^0, r) 

|L >P ) = /|t_, P )- fl |t+, P ) 



(7a) 
(7b) 
(7c) 
(7d) 



in real space with u = cos(#), v = sin(0), / = cos(<^) 
and g = sin(<^), the position r = (r, s) and the wave 
vector Qo = (tt, 0) of the minimum of the dispersion. 
The triplon states \t m ) with m G { — ,0,+} are given 
by \t-) = l/y/2(\t x )-i\t y )), \t ) = \t z ) and \t+) = 
\/y/2(\t x ) + i \ty)). The tensor product of all singlet states 
\s r ) is the vacuum |0), so that the hardcore triplon creation 
operator with m G { — , 0, +} is defined by r |0) := |t m>r ) 
and the annihilation by t m ^ r |t m>r ) := |0) and so on. In 
this basis the magnetic field is described by the operator 
— h(t\t + —t[_t_). The two independent variables and 
if are varied to minimize the classical ground state en- 
ergy. This choice also ensures that (i) all linear terms in 
the triplon operators vanish and (ii) a massless Goldstone 
mode appears as it has to be. 

Previous work [38. 39 ; 51 applied the transformation ([7]) 
to the original spin model. This is not possible for IPA- 
CUCI3 because the dimers are too strongly coupled. Hence 
the CUT is mandatory and we apply the real space trans- 
formation (|7]) to #iD,eff + Winter, eff from Eqs. (|2|4p keeping 
the bilinear terms. Fourier transformation and Bogoliubov 
diagonalization finally provides the lowest lying modes. 
Their resulting gap energies are displayed in Fig. HJ No 
parameters are adjusted. 
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Fig. 4: (Color online) Gaps in IPA-CUCI3 vs. the reduced mag- 
netic field gH/2. Solid lines show theoretical results, see main 
text; symbols mark experimental data from Refs. [23] (setup I 
& IV) and 9^ (setup II & III). 

Second, the upper critical field H C 2 can be determined 
exactly for the spin model ([Taj) to be H C 2 = (2J2 + 
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J$)I{9V>b) ~ 2/flf • 45.8 T. After the transformation © 
is applied to the dispersion obtained from CUT we obtain 
H C 2 « 2/g - 45.1 T. The very good agreement of these 
two values strongly supports the approximations made. 
Additionally, the theoretical values also match the exper- 
imental result [31] H c2 = (43.9 ± 0.1) T(2/g) within 4%. 
In view of the neglect of anisotropics and magnetoelastic 
effects, cf. Ref. [22], this nice agreement lends independent 
support to the advocated microscopic model. 



as done previously [40, 53 because IPA-CUCI3 is not far 
enough in the dimer limit. 

For comparison, we also include A(T) as derived from 
the nonlinear a model on 1-loop level [55] in Fig. [5] (dashed 
line). It is obtained from 
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Fig. 5: (Color online) Spin gap in IPA-C11CI3 vs. temperature 
T. Lines show theoretical results from CUTs and mean-field 
(solid) and from the nonlinear a model (dashed), respectively. 
Inset: Temperature dependence of the condensate fraction s 2 . 

Third, the temperature dependence of the gap A(T) 
supports that the low-lying excitations are hardcore 
triplons. We apply the mean- field approach in Refs. 
P0B6ll52H 5l] to #iD,eff + Winter, eff from Eqs. In each 

nonlocal term (4, r t m?r , or 4,r*-m, r ' or W*-m,r' with 
r 7^ r') all creation operators r are multiplied by the 
singlet annihilation s r and the annihilation operators t m>r 
by the singlet creation sj. Local terms remain unchanged 
because they do not change the local singlet number. Fi- 
nally all singlet operators are replaced by the condensate 
value s(T) = (s^) = (s) with s G [0,1]. In a nutshell, a 
factor s 2 appears in front of each nonlocal term. 

This implies a dependence of the dispersion on s and 
hence on temperature [IU][52j[5l], denoted by c<j s ( T ) (ft, /). 
The self-consistent solution is found from the hardcore 



condition 1 = (sjs r - 
1 - 32/(1 + 32) with 



EmiAr) leading to s 2 (T) = 



rl/2 



pl/Z nL/ 

/ dft / dl e -P"s ( T)(Ki)^ ^ 
J -1/2 J -1/2 



Figure [5] compares the result (solid line) of this simple 
approximation to INS data [32|l33] . Up to 15 K the exper- 
imental data is matched perfectly. We attribute the dis- 
crepancy at higher temperatures to the insufficient treat- 
ment of the hardcore constraint by the above approach 
(for 15 K the condensate fraction s 2 is only 0.77). Note 
that we only apply the mean-field theory to the disper- 
sion obtained from CUT, not to the original spin model 



with u(h,l,T) := y/u 2 (h, I) + A 2 (T) - A 2 (0); the con- 
stant C is determined for T = 0. Interestingly, this ap- 
proach describes the experimental data less accurately if 
the experimental dispersion at T = is used for a; (ft, Z), 
cf. Ref. [33]. We presume that the hardcore constraint is 
not accounted for sufficiently well by Eq. ©. 

In summary, we showed that the available experimental 
evidence for IPA-CUCI3 is consistent with a quantitative 
model of weakly coupled asymmetric S = 1/2 spin ladders 
with hardcore triplons as excitations. Such systems are of 
great current interest because they allow for the study of 
Bose-Einstein condensation of triplons and of the mass- 
less excitations above this condensate [8,9.. Additionally, 
they represent gapped quantum liquids known to display 
considerable quasiparticle decay [T3HT6] . 

Our high-precision analyses of inelastic neutron scatter- 
ing data and of the temperature dependence of the mag- 
netic susceptibility is based on advances in continuous uni- 
tary transformations [16] and high temperature series ex- 
pansions. The established quantitative model paves the 
way for further quantitative studies, both experimental 
and theoretical, of the decay of massive quasiparticles and 
of the condensation of hardcore bosons. 

The latter is illustrated by the excellent agreement of 
the calculated gap energies as function of magnetic field. 
Additionally, the description of IPA-CUCI3 by dispersive 
hardcore triplons is strongly supported by the agreement 
of the temperature dependence of the spin gap. 

By this work, a quantitative model for IPA-CUCI3 is 
established. Concomitantly, we exemplarily showed how 
CUT results in one dimension at zero temperature and 
zero magnetic field can be extended to render a quanti- 
tative description in two dimensions at finite temperature 
and finite magnetic field possible. We expect this approach 
to continue to be fruitful also for other systems. 
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